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Abstract 

Theoretical calculations in the framework of Kane model have been carried 

out in order to elucidate the role of interband mixing in forming the energy 

spectrum of two-dimensional carriers, localized in a surface quantum well in 

narrow gap semiconductor. Of interest was the mixing between the 2D states 

and heavy hole states in the volume of semiconductor. It has been shown that 

the interband mixing results in two effects: the broadening of 2D energy levels 

and their shift, which are mostly pronounced for semiconductors with high 

doping level. The interband mixing has been found to influence mostly the 

effective mass of 2D carriers for large their concentration, whereas it slightly 

changes the subband distribution in a wide concentration range. 
PACS number(s): 73.20.At, 73.40.Qv, 73.40.Gk 
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I. INTRODUCTION 



Interband mixing effects attract much interest because they can play a decisive part in 
forming the energy spectrum of low- dimensional systems. Such is in the case for InAs-GaSb 
semimetallic superlattices, whose physical properties are determined by the mixing between 
valence band states of GaSb and conduction band states of InAs at the interfaces. |T]-§[ 

Analogous situation takes place in a metal-insulator-semiconductor structure based on a 
narrow gap semiconductor. A small value of the energy gap in the semiconductor leads 
to the fact that at energies below the top of the valence band the two-dimensional (2D) 
states localized in a surface quantum well can mix with the valence band states and 2D 
carriers can tunnel from the space charge region into the valence band states in the volume 
of semiconductor. This results in a shift of the 2D subband energies and a broadening of 
these levels. 0-0 Due to high value of the effective mass, the probability of tunneling of 2D 
carriers into the heavy hole states is small. Therefore the interband mixing with the heavy 
hole states was neglected. This approach is justified for semiconductors with low impurity 
density, because the surface potential in such materials is smooth in the shape, and the 2D 
states localized in the quantum well are spatially separated from the valence band by a wide 
region of forbidden energies. 

In semiconductors with high doping level the effects of interband tunneling of 2D electrons 
into the heavy hole states can be stronger. The influence of these effects on the energy 
spectrum and broadening of the 2D states localized in the surface quantum well in the 
narrow gap semiconductor (HgCd)Te is discussed in the present article. 

II. THEORETICAL MODEL 

A small value of the energy gap (Eg) in semiconductor investigated makes it necessary 
to employ a multiband Hamiltonian in the fcP-based calculations of the energy spectrum. 
We started from the multiband Hamiltonian derived from Kane's model making the usual 
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assumption that the energy difference between the valence Fg-band and spin-orbit spht-off 
Fy-band is infinite. To consider the effects involving the heavy hole states, the interaction 
with remote bands has been taken into account by a standard procedure of including ad- 



ditional terms with 7-parameters in the Hamiltonian. [|T0[ An isotropic approximation has 
been used. 

The Kane Hamiltonian, which is a 6 x 6 matrix in the above assumptions, can be block 
diagonalized, when the axis z is chosen to be normal to the interface and the axis y along 
the direction of the carrier motion. In this case the Hamiltonian is two 3x3 matrices for 
two groups of states. One of the matrices is defined by 



3 \ 2 dz J 



k a 

'2 dz 



g = E^Hz) + 



2mn 



(71 + 27) 



dz^ 



(71-7)^^ +e(fi{z), 



T = E'Hz) + 
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(71 - 27) 



(71 + 7)^^ + eip{z), 



where k = ky, P is the momentum matrix element, 71, 7 are the parameters, which 
describe the interaction with remote bands, ^p{z) is the electrostatic potential, and E'^'^{z), 
E^^{z) are the energies of the conduction and valence band edges, respectively. Here, the 
assumption has been made that the values of P, 71, and 7 are independent of coordinate. 

The second matrix H~ is obtained from by replacing k by —k. 

The electrostatic potential can be derived from the Poisson equation. We assumed ^p{z) 
to be parabolic in the space charge region and constant in the volume of semiconductor: 



^s{l - z/Lf, 0<z<L 



0, 



z > L, 



(2) 



where 
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Here, Na — Nd, and n stand for density of acceptors uncompensated, and dielectric constant, 
respectively. Expression is a good approximation for the case when the concentration of 
2D electrons is less than the value of {N^ — Nij)L. 

Thus, Schrodinger equation is two independent systems of differential equations of the 
second order. For solving these systems boundary conditions have to be deduced. The 
introduction of the boundary conditions is based on a model proposed in Refs. |Q and fill 



for calculation of the energy spectrum of 2D states in MIS structures, based on ordinary 
and inverted semiconductors, respectively. As in these papers, we supposed (see Fig.|^a): (i) 
the band structure of insulator to be identical to that of semiconductor; (ii) the energy gap 
in the insulator to be much greater than that in semiconductor, and the conduction and 
valence band offsets, denoted as Dc and D^, respectively, to exceed the value of Eg] (iii) the 
energy bands in the insulator to be "flat", i.e. (f{z < —d) = const, and (iv) the structure to 
have a smooth insulator /semi conductor interface of width d with linear dependence of the 
band-edge energies on z coordinate. |TT| 

Since a general solution of the Kane Hamiltonian for (p{z) = const is known, the wave 
function on the insulator {z < —d) and semiconductor {z > L) sides of the structure can be 
simply derived. 

For z < —d it contains only the terms exponentially vanishing deep into the insulator. 
To choose the boundary conditions on the semiconductor side let us briefly consider three 
different cases, which can be realized depending on the energy and longitudinal component 
of quasimomentum of 2D state. First, the 2D state is not degenerate with the valence band 
(region I in Fig.^D). In this case the wave function exponentially diminishes deep into the 
volume of semiconductor. Second, the 2D state is degenerate only with the heavy hole 
valence band (region II). In this case the normal component of the light hole quasimomen- 
tum is imaginary. Therefore, the wave function in this case is a superposition of the light 
hole term, which diminishes exponentially deep into the semiconductor, and two oscillating 
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terms presenting the plane wave associated with the heavy hole. Third, the 2D state is in 
resonance with both heavy and light hole valence bands (region III) and the wave function 
is a superposition of two plane waves corresponding to the light and heavy holes. In the 
present paper we focus our attention on the second case. 

Thus, starting from the exact solution of the Kane Hamiltonian on the insulator side 
we numerically integrated the differential equation systems through the space charge region 
and chose those solutions which satisfy the boundary conditions on the semiconductor side. 

It is clear that in region II the eigenvalue problem has a continuous spectrum of solutions, 
i.e. at any energy value we can find the wave function satisfying the boundary conditions. To 
define the energy level associated with the resonant 2D states and its broadening caused by 
the interband mixing between the 2D states and heavy hole states, the Levinson's theorem 
has been applied. 



III. DISCUSSION 

To analyze peculiarities of the energy spectrum of 2D states arising from the interband 
mixing, we have performed numerical calculations using realistic parameters which are close 
to the parameters of the sample investigated in Ref. [Q. They are the following: Eg = 50 
meV, A^^ — = 6 x 10^^ cm~'^, k = 20. The calculations have been carried out using 
two sets of 7-parameters: 71 = 2.0, 7 = (Ref. [0]) and 71 = 4.5, 7 = 1.0 (Ref. 
To make identical the dispersion law of bulk electron states calculated with different sets 
of 7-parameters, we changed slightly the value of the momentum matrix element P from 
8.0 X 10~^ eV-cm, for the first set, to 8.1 x lO"*^ eV ■cm, for the second one. Parameters of 



the insulator/semiconductor interface are the same as in Ref. Dc = 2 eV, Dy = 1 eV, 
and d = 10 A. 

Figure ^ shows the energy spectrum of 2D states calculated with different sets of 7- 
parameters as compared with the energy spectrum of 2D states calculated in the infinite 
heavy hole mass approximation with a so-called "mid-gap" boundary condition on the semi- 
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conductor side. [|T3| Two branches of the energy spectrum specified by and k~ are the 
results of calculation for two groups of states. These branches are two "spin" branches of the 
ground 2D subband split by spin-orbit interaction in an asymmetric quantum well. |I^-[T8[] 
It is clearly seen from the figure that the interband mixing of 2D states with the heavy hole 
states results in the energy shift of 2D subbands. The smaller is the heavy hole effective 
mass {rrih = "^0/(71 — 27)) the lower is the energy position of the 2D subbands. 

One more important feature of the energy spectrum of 2D states appearing due to in- 
terband mixing is the broadening of 2D energy levels. This broadening is shown in Fig.|^ as 
hatched region. As is seen the broadening of 2D levels is negligibly small for energies close 
to the top of the valence band in the volume of semiconductor. At more negative energies 
both 2D "spin" sublevels become broad. The value of broadening is different for and k~ 
sublevels for a fixed energy. It increases with decreasing energy or quasimomentum value. 
The maximum value of the broadening is about 10 meV. It is comparable to the broadening 
of 2D states appearing due to tunneling into the light hole states which has been calculated 
in Ref. ^ for a sample with close parameters. On further decreasing quasimomentum the 
2D states fall in resonance with the light hole valence band (region III in Fig. [^). Note 
that the broadening value slightly depends on the heavy hole effective mass. 

Thus, the interband mixing between the 2D states and heavy hole states leads to two 
effects: the shift of subband energies and the broadening of 2D levels. Let us consider how 
these effects can manifest themselves in experiment. 

Traditional experimental methods, such as galvanomagnetic investigations, volt- capac- 
itance spectroscopy, give information about the carriers at the Fermi energy. In p-type 
(IlgCd)Te the Fermi level lies near the top of the valence band in the energy range —5. ..+5 
meV at low temperatures. As was shown above for these energies the broadening of 2D 
states is small and therefore these methods must not be sensitive to this effect. Tunneling 
investigations allow to get information about the 2D states in a wide energy range including 
the energies, where the broadening is large. Besides, the maximum value of broadening is 
comparable with the cyclotron energy in magnetic fields right up to 2 T. This can result 
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in additional broadening of resonance peaks in magneto-optical and cyclotron resonance 
experiments. 

As to the effect connected with the influence of interband mixing on the dispersion law of 
2D states traditional experimental methods give no way of direct deducing a quasimomen- 
tum dependence of the energy of 2D states for a fixed value of a surface quantum well. For 
example, in tunneling experiments [0,0 a voltage bias, applied to the structure, changes 
not only the energy tested, but surface potential as well. The value of the surface potential 
is usually used as a fitting parameter in analyzing experimental data, for example, depen- 
dences of 2D subband concentration and the effective mass of 2D carriers on their total 
concentration. Figure ^ shows the relationship between 2D subband concentration and the 
total concentration of 2D electrons, calculated using two sets of 7-parameters and infinite 
heavy hole effective mass. The concentration of 2D electrons is changed by varying the 
surface quantum well depth from 250 to 320 mV. It is clearly seen from the figure that 
the interband mixing does not cause essential changes in the subband distribution of 2D 
electrons. One can see only slight increase in ratio between concentrations of 2D electrons 
in and k~ "spin" subbands with decreasing the heavy hole effective mass. Thus, in spite 
of the fact that the interband mixing results in shifting the energy positions of 2D subbands, 
it does not strongly influence the subband distribution of 2D electrons. 

The effective mass of carriers is a differential characteristic of the energy spectrum and 
it might be more sensitive to the interband mixing effects. We have calculated the effective 
mass of 2D electrons m* = h^k {dE/dk)~^ at E = Ep function of their total concentra- 
tion (Fig. I). As is clearly seen for large concentration of 2D carriers (A^^ > 0.55 x lO^^cm"^) 
the effective mass of 2D electrons, calculated in the model taking into account the interband 
mixing (solid and dashed curves), is larger than that calculated with the infinite heavy hole 
mass (dotted curves). The difference increases with increasing the total concentration of 2D 
electrons. It is evident that the concentration dependence of the effective mass, calculated 
with 7i = 4.5, 7 = 1.0 deviates mostly from the results obtained in the infinite heavy hole 
mass approximation. For /c+-"spin" states this difference is larger and at Ns — lO^^cm"^ it 

7 



comes to more than 10%. 

A distinctive situation occurs for small concentrations of 2D electrons (A^^ < 0.45 x 
lO^^cm^^) for fc^-group of states. The effective mass, calculated with 7-parameters, mono- 
tonically decreases with decreasing the total concentration, whereas the effective mass, cal- 
culated in the framework of the infinite heavy hole mass approximation, goes through a 
minimum, and at some value of the total concentration it becomes larger than the effective 
mass calculated with 7-parameters. To explain this peculiarity, let us turn to Figure |^. One 
can see that the bottom of fc^-subband lies at 7^ (the dotted curve in the figure). As a 
consequence, the effective mass m* = ti^k {dE/dk)^^ turns to infinity on the bottom and is 
sufficiently large in its vicinity. For the low concentrations of 2D electrons, this portion of 
the energy spectrum is found to be near the Fermi level for the case of infinite heavy hole 
mass, whereas it is lower in the energy when the finite value of the heavy hole mass is taken 
into account. This reflects in the behavior of the effective mass of 2D carriers for their low 
concentration. 

Thus, the presented results show that the interband mixing causes the effective mass of 
2D states to be larger than that calculated in the infinite heavy hole mass approximation. 
This is more evident for fc"'"-"spin" states and for large concentration of 2D carriers. The 



effective mass of 2D carriers calculated with the most-used set of 7-parameters ||12| differs 
mostly from the results of calculation, in which the interaction with remote bands is ne- 
glected. To our knowledge there are no experimental data concerning the effective mass 
of 2D states in highly doped narrow gap semiconductors in case, when 2D states are in 
resonance with the heavy hole valence band. 

IV. CONCLUSION 

The effects of interband mixing between 2D states and heavy hole valence band states 
are found to influence the energy spectrum of 2D states localized in a surface quantum 
well in the narrow gap semiconductor with high doping level. They manifest themselves in 
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the broadening of 2D energy levels, which is sufficiently large for the energies close to the 
bottom of 2D subband, and the shift of 2D subbands to lower energies. The latter does not 
practically affect the concentration dependence of subband distribution of 2D carriers, but 
leads to a change of their effective mass for high total concentrations of 2D carriers. 
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FIGURES 

FIG. 1. (a) Model of an insulator - narrow gap semiconductor structure with a surface quantum 
well used in the calculations, (b) The dispersion law E{k, kz = 0) for carriers in the volume of 
semiconductors (solid curves). The hatched regions show the continuous spectrum. The numbers 
I, II, and III denote three different regions of the energy spectrum of 2D states (for details, see 
text) 

FIG. 2. The energy spectrum of 2D states calculated with different sets of 7-parameters: 
71 = 4.5, 7 = 1.0 (the solid curves) and 71 = 2.0, 7 = (the dashed curves), as compared 
with the energy spectrum of 2D states calculated in the infinite heavy hole mass approximation 
(the dotted curves) . The value of the surface potential used in the calculations is 285 mV. A;+ and 
k~ denote two different "spin" branches of the energy spectrum. The dot-dash curves show the 
borders of region II, where the 2D states are in resonance with the heavy hole valence band. 

FIG. 3. The dependence of the electron concentration in each 2D subband on the total concen- 
tration of 2D electrons. The notations and curve types are the same as in Fig. 2. 

FIG. 4. The dependence of the electron effective mass on the total concentration of 2D electrons. 

The notations and curve types are the same as in Fig. 2. 
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